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Abstract — Uncertainties in the retrievals of microwave 
I a nd -surface emfssivities are quantified over two types of 
land surfaces: desert and tropical rain forest. Retrievals from 
satellite- based microwave Imagers, including the Special Sensor 
Microwave Imager, the Tropical Rainfall Measuring Mission 
Microwave Imager, and the Advanced Microwave Scanning 
Radiometer for Earth Observing System, are studied. Our 
results show that there are considerable differences between the 
retrievals from different sensors and from different groups over 
these two land-surface types* In addition, the mean emissivity 
values show different spectral behavior across the frequencies* 
With the true emissivity assumed largely constant over both 
of the two sites throughout the study period, the differences 
are largely attributed to the systematic and random errors in 
the retrievals. Generally, these retrievals fend to agree better 
at lower frequencies than at higher ones, with systematic 
differences ranging l%-4% (3—12 K) over desert and l%-7% 
(3-20 K) over rainforest. The random errors within each 
retrieval dataset are in the range of 0.5% -2% (2-6 K). In 
particular, at 85 *5/89.0 CHz, there are very large differences 
between the different retrieval datasets, and within each retrieval 
dataset itself. Further investigation reveals that these differences 
are most likely caused by rain/cloud contamination, which can 
lead to random errors up to I CM 7 K under the most severe 
conditions* 

Index Terms — Brightness temperature, land-surface emissivity, 
measurement uncertainty, microwave radiometry, random errors, 
remote sensing, systematic errors. 
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L Introduction 

L AND surface emissivity at microwave frequencies eon- 
tains a wealth of information on the physical, biological, 
and hydrological stales and processes of the earth’s surface. 
This forms the basis lor remote sensing of a wide range of 
land-surface states and processes such as soil moisture, vegeta- 
tion characteristics, and land-cover dynamics 1 1 ]. In addition, 
land-surface emissivity acts as the background signal for the 
retrieval of atmospheric variables, such as water vapor, rainfall, 
and snowfall, and therefore greatly affects the accuracy and 
uncertainty in such measurements [2], [3]. 

Satellite-based emissivity retrievals have been performed for 
nearly three decades [4]-[9], while the methodology remains 
largely unchanged. Microwave emissivilies arc derived from 
satellite- based observations though radiative transfer tabula- 
tions* Microwave radiometers onboard polar-orbiting satellites 
produce brightness temperature (Tb) measurements at the 
top of the atmosphere (TO A). With atmospheric tempera- 
ture and moisture profile data, one can remove the portion 
of Tb originated from the atmosphere, to obtain the microwave 
emission from the land surface. Subsequently the land-surface 
emissivity can be computed if the surface temperature effective 
for the emission is known. More recently, the retrieval of 
emissivity has also been implemented in a variational [10] 
and/or iterative [11] framework in which many variables 
affecting the radiative transfer processes, including emissivity, 
can be estimated simultaneously* 

Usually, emissivity retrievals are only performed for clear 
days, due to the difficulty in estimating the atmospheric 
contribution from a cloudy or rainy atmosphere, and to the 
strong atmospheric scattering and absorption of land-surface 
signals under such conditions, especially at higher frequencies. 
However, even for a cloud-free atmosphere, there are many 
error sources that lead to uncertainties in emissivity retrievals, 
including instrumental errors, inaccuracies in the atmospheric 
profile data, imperfect cloud screening, and misrepresentation 
or the land-surface temperature [7J, [14], 115]* In addition, the 
heterogeneity of the land -surface radiometric properties and 
the shills in instrument footprint locations introduce sampling 
errors that enhance the uncertainties. 

Despite its importance, the uncertainty in emissivity 
retrievals has not been well quantified. This is reflected in 
the considerable disagreements among the most recent, state- 
of-the-art datasets |16|. The leading difficulty is the lack of 
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TABLE 1 

FREQUENCIES OF THE MICROWAVE IMAGERS 


frequencies of Microwave Imagers (GHz) 


SSM/l 



19.35 

22.2 (V} 

37.0 

85.5 

TMI 


1 0.65 

19.35 

21.3 (V) 

37.0 

85.5 

AMSR-E 

6.9 

m.65 

lft.7 

2.3.S 

36.5 

89.0 


“ground-truth” data, especially on Ihe global scale, Most of 
the Held campaigns for land emissivity studies are short-lived 
and small-scale ones, and generally they are not carried out 
in coordination with any specific satellite-based instruments 
or overpasses. Without reliable reference data, the immediate 
impediment to uncertainty quantification is the inability to 
apportion the variability to measurement error or to natural 
variability of the land-surface emissivity. 

In this paper, we circumvent this difficulty by strategi- 
cally selecting two types of land surfaces whose emissivities 
are largely constant: the Sahara Desert and the Amazon 
Rainforest. Then all the variations within a set of retrievals 
are caused by the uncertainties, defined as the spread among 
independent measurements, Although this paper is limited to 
these two types of land surfaces, it represents a practical effort 
to make progress in solving an otherwise intractable problem. 
Thus, this paper not only provides the first quantitative results, 
but also facilitates more educated inference on the magnitude 
of uncertainty over other surfaces. 

The emissivity datasets and methodology we employed are 
described in the following section. In Section III, we provide 
observational evidence to substantiate our assumption of the 
constant emissivity over the desert and rainforest areas. Rased 
on such an assumption, we present results in Section IV to 
quantify the uncertainties. The results are then summarized 
and discussed in Section V. 

II. Data and Methodology 

To support the upcoming Global Precipitation Measurement 
(GPM) mission, NASA’s Precipitation Measurement Missions 
(PMM) Science Team formed the Land-Surface Working 
Group (LSWG) to improve land-surface characterization at 
microwave frequencies. LSWG has assembled a collection of 
clear-sky land-surface emissivity retrievals from many contem- 
porary space- borne passive microwave sensors, over selected, 
representative land-surface types such as desert, rainforest, 
midlatitude agricultural land, wet land, and high-latitude cold 
regions |16|. The collection of sensors includes the vSpecial 
Sensor Microwave Imager (SSM/1), the fro pi cal Rainfall Mea- 
suring Mission (TRMM) Microwave Imager (TMI), WindSat 
aboard the Coriolis satellite, the Advanced Microwave Sound- 
ing Unit (AMSU), and the Advanced Microwave Scanning 
Radiometer for Earth Observing System (AMSR-E), This data 
collection greatly facilitates intcrcomparison and evaluation of 
land- surface microwave retrievals, and enables us to assess the 
current skills and uncertainties. 

In this paper, we focus our evaluation on the three conical- 
scan microwave imagers: SSM/1, TMI, and AMSR-E (Table I), 
because they have a similarly wide frequency range, and most 
of their channels have both vertical and horizontal linear polar- 
izations. Instantaneous retrievals, including both ascending and 


descending passes, are used whenever available. Retrievals 
from different data providers for the same imager are all 
included as independent estimates, as they are mostly derived 
from different algorithms and/or with different ancillary data. 
A common one- year period — July 1, 2006 to June 30, 
2007 — is used for our paper. 

For SSM/I, retrievals from three Defense Meteorological 
Satellite Program (DMSP) platforms (FI 3, FI 4, and FI 5) arc 
included. The data are provided by the Centre National de 
la Recherche Sctentifique (CNRS), France, and the retrieval 
procedure is described in [6J and |8|. TMI retrievals are 
produced by Nagoya University, Japan, with a similar method. 
For AMSR-E, two independent retrieval datasets are used. 
One is produced by NOAA’s Microwave Integrated Retrieval 
System (M1RS) using a 1-D variational algorithm [10], and 
is denoted as AMSR-E (MIRS). The other is produced by 
NOAA’s Cooperative Remote Sensing Science and Technology 
Center (CREST), and the retrieval procedure is documented in 
[17] and [18J. This dataset is denoted as AMSR-E (CREST). 

Uncertainty gauges our ignorance and state of knowledge* 
In practice, uncertainty can be quantified as the spread or 
disagreement among independent measurements of the same 
physical quantity. If each of the measurements has no system- 
atic biases, then such spread arises solely from the random 
errors. However, most often, each measurement has distinct 
systematic biases that contribute to the total uncertainty. 
Difficulty arises when there are no “ground -truth” data avail- 
able: one is not able to separate systematic errors from random 
ones, and cither type of error from the natural variability of 
the measured quantity. Under these conditions, uncertainly 
quanti lication is impossible, 

A practical starting point lo elude this difficulty is to 
quantify the uncertainty over areas where the physical 
variable — here microwave emissivity — js constant. This is the 
one case in which “ground- truth” data are not needed, as all 
the variations in the data are from the measurement errors. 
One may still not be able to identify the absolute amplitude 
of systematic errors, but their differences among independent 
measurements can be obtained, and these differences provide 
substantial insight into their reliability as an ensemble. More- 
over, the random errors can be easily quantified, which more 
often arc the dominant pari of the total uncertainty. 

For this paper, we selected two land-surface types, the 
Sahara Desert and the Amazon Rainforest. In these regions, 
the desert roughness and the vegetation properties of tropical 
forest are not expected to change significantly with time. 
Since roughness and vegetation are two key controlling factors 
for microwave emissivity at these frequencies, they can be 
used as a constant- reference surface. This property has been 
exploited for validation and calibration of passive microwave 
sensors and other land-surface parameter retrievals [19J— [27J. 
In the following section, we will substantiate this assump- 
tion with the long-term AMSR-E observations. Though the 
true value of the assumed constant emissivity is unknown, 
we can compare the independent retrieval datasets and their 
disagreements to infer the magnitude of the uncertainties, 
including both systematic and random errors. The two sites 
we used in this paper are designated by LSWG as “Desert” 
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and “Amazon!” located at (22 fl N, 29° Ll) and (2° N, 55° W), 
respectively. 


and since ' \ C l \ 2 for most land surfaces, the above 
equation can be approximately rewritten as 


111, Constant Emissivity Over Study Sites 


To prove that the land-surface emissivitics are approx- 
imately con slant over the Desert and Amazon2 sites, we 
examined the microwave polarization difference index (MPDI) 
over these two Locations. MPDI is essentially a normalized 
measure of the polarization. Ilie MPDI is computed from 
the daily AMSR-F Ih data over the three-year period of 
July 2004 through June 2007, for all weather conditions. 
Tb-based MPDI is defined as [28]-[31] 


MPDli D 


Tbv Tbh 
Tb v C Tbh 


where Tbv and Tbh are the TOA Tb of vertical and horizontal 
polarization, respectively, for a specific frequency and scan 
angle. Similarly, emissivity -based MPDI is defined as [13J 

■t u 

MPDI" D 

v ^ h 


where " v and "j, are the land-surface emissivity of vertical and 
horizontal polarizations, respectively, for a specific frequency 
and scan angle, 

MPDI is closely related to the land-surface states. 'Hits 
is because with u dear-sky atmosphere, the polarization 
difference originates exclusively from the land surface; the 
atmosphere's attenuation and emission generate no additional 
polarization, serving only to suppress the polarized signals 
down below. Thus, Tb-bascd MPDI, which contains both 
signals from the fund surface and (unpolarized) signals from 
the atmosphere, cun be shown to be a lower-bound estimate 
of the emissivity-based MPDI. 

For a nonsc altering, plane- parallel atmosphere and for a 
given scan angle and frequency, ihe TOA brightness temper- 
atures for vertical and horizontal polarizations can be com- 
puted by the following integrated radiative transfer equations 
[12], [131V 

Tb v D T,“ C T*. I U " v /e ac C V' v e uc 

D T" C Tjfe 1 C ,T S C T a # /V nn 0) 
Tbh D T™ C T*. I O "h/e DC C T t \e ac 

D T» C Tfe 1 C .T s G T*l\e JC (2) 


where T a and T, are the upward and downward (including 
cosmic background) radiations through the atmosphere, T, the 
land-surface temperature, and the atmospheric optical depth 
along the view path, respectively. Then from (1) and (2) one 
can derive 


MPDI t D 
D 
D 


Tby Tbh 
Tb v C Tb h 

,T t DT”/e r ," v ay 

2.Tj C Tj # e 1C .% n T a # /e‘ C \l 
MPDI- 

I q 2.T"CT*f I 
T,! 7, It \C\! 


(3) 


mpdi t 


MPDI- 


I C 


T a ~ e CTjf 
T, T* 


(4> 


As T 5 > T a for Ihe earth's clear atmosphere in the fre- 
quency range under siudy, we have MPDIt - MPDT . In other 
words, the Tb-bascd MPDI is a lower-bound estimate of the 
emissivity -based MPDI. The more transparent the atmosphere 
is (the smaller the values of T a ,T a # , and are), the closer the 
two MPDIs will become | I3J. 

Fig. 1 shows the Tb-hased MPDI for all the six AMSR-E 
channels, over the Desert [Fig. 1(a)], Amazon2 [Fig. 1(b)] and 
a third LSWG site, the southern great plains [SGP; Fig. 1(c)]. 
Over l he three-year period, the MPDI for both Desert and 
Amazon! remains fairly constant. In contrast, SGP’s MPDI 
values exhibit strong seasonal variations across all the channels 
that are related to the seasonal variations in vegetation and soil 
moisture. Though the data contain both clear and cloudy days, 
the latter only serve to further reduce MPDI values and do 
not obscure variations during clear days. Comparing Fig, t (a) 
and (b) with (c), we can verify that our assumption of constant 
emissivity over Desert and Amazon2 is valid. 

Over desert, the MPDI values are large, and so are the dif- 
ferences among the channels. This is consistent with the wciS- 
knowm strong polarization signature over desert areas [32], 
Nevertheless, over the three-year period, the values are remark- 
ably constant, especially for lower frequencies, which are 
not sensitive to the atmosphere. Slight seasonal variations are 
present in the water vapor channel (23.8 GHz), which suggests 
they are mostly atmosphere- induced. The synchronized, less- 
pronounced fluctuations in the 89.0-GHz channel indicate their 
atmospheric origin as well. 

Over Amazon2 | Fig. 1(b)], the MPDI values are much 
lower, and close to zero for all the frequencies. This is con- 
sistent with the fuel that over ihe dense forest the microwave 
signals are not. strongly polarized. Despite the dense canopy, 
the lower frequencies exhibit slightly higher MPDI values, 
indicating some of the emission from either exposed soil 
surface or water bodies. Indeed, satellite images from Google 
Barth ( h ttp ://earth . google . com ) indicate there is a river nearby 
well within the sensors 1 footprints, Ihe higher MPDI values 
are also consistent with results shown in Njoku et al.[24|. 
Overall, one can see MPDI remains fairly constant and 
very small across all the channels. Thus, it is reasonable to 
assume the emissivity values over Amazon2 are approximately 
constant loo. 

Wc also examined another arbitrarily selected, distant desert 
site (22° N, 5 W W) and Amazon site ( T S, 70° W). The 
temporal variation in their respective MPDI time series is 
similarly small (not shown), indicating the slationarily of 
MPDI over these two types of surfaces a robust feature. We 
expect that the other imagers (SSM/I and I Ml) would show the 
same behavior, based on results from intersensor comparison 
studies [l9]-[27]. 
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rig. i. MPDI calculated from AMSR-E brightness Lei literatures (Tb) fur ascending passes over a three-year periud (July 2004 through June 2007), over the 
three LSWG sites (a) Desert, (b) Amazon2, and (c) SGP. The cause of the discontinuities in (a) near October 2004 and April 2005 is not documented but 
they are confined to this particular region. 


IV. RESULTS 

In this section, we first present Itie basic characteristics 
of both the systematic and random errors, including their 
statistical distributions shown as histograms. To gain additional 
insight, we then examine the errors in MPDI space, with 
comparisons with AMSR-B IT} data. Finally, we provide 
diagnostic analysis on the sources of some of the error features 
identified in this paper. 

A Systematic and Random Errors 

Land- surface emissivity retrievals over the two evalua- 
tion sites show considerable systematic and random errors. 


Fig. 2 shows box -and- whisker plots for the six datasets 
and for both polarizations. Over the Desert site [Fig. 2(a) 
and (b)], the ensemble as a whole showed the expected 
behavior — vertically polarized emissivity decreases with fre- 
quency [Fig. 2(a)], while horizontally polarized emissivity 
increases [Fig, 2(b) J [32 J. However, there are considerable 
systematic differences between the mean values from each 
retrieval dataset, indicating most of them, if not all, have 
systematic errors, regardless of the unknown "ground truth.” 
T he systematic differences are the smallest at lower frequen- 
cies (6.9 and 10.65 GHz), reflecting partly their insensi- 
tivity to atmospheric effects. T he random errors also show 
a strong dependency on frequency. T he higher frequencies 
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Fig. 2. Intercomparison of microwave emissivity retrievals at two LSWG sites Desert (lop) and Ama/.on2 (bottom), over a one-year period from July 1, 
2006 to June 30, 2007, for SSM/T (FI 3, FI 4, and FI 5), TM1 and AMSR-E. Both vertical (left) and horizontal (right) polarizations are shown. The extremely 
high values in the vertically polarized channel of SSM/1 F 1 5 at 22.2 GHz. are caused by an instrument problem documented in [33]. (a) Desen, V-pol. 
(b) Desert, H-pol. (c) Amazon2, V-pol. (d) Amazon2, H-pol. 


(85.5/89.0 GHz) tend to have the highest spread from their 
mean values, suggesting again atmospheric effects are the 
source. The vertically polarized channel of SSM/I FI 5 at 
22.2 GHz shows extremely high values [Fig. 2(a) and 
(c)l This is caused by an instrument problem documented 
in [33]. 

Over Amazon2 | Fig. 2{a) and (d)j, the magnitude of the 
systematic differences among the retrievals is similar to that 
of the Desert site, except that the differences at the higher 
frequencies (85.5/89.0 GHz) are much larger. Similarly, the 
random errors are also much higher at these frequencies. I bis 
also suggests that the atmospheric effects are playing an even 
larger role here, considering that Amazon2 has a much moister 
atmosphere and many more cloudy/rainy days than those of the 
desert site. 

There is a lack of smoothness in the mean emissivity 
spectra at either site. For example, CREST’s AMSR-E 
retrievals fashioned a bump at 23,8 GHz in its horizontal 
polarization at the Desert site, but it has a dip at 36-5 
GHz in both polarizations at Amazon2. Other retrievals 
show such bumpiness in varying degrees. We believe the 
emissivity spectra should be smooth and monotonic, because 
over the frequency range being studied the land surface 
does not have any known physical mechanism that responds 
differently to a particular frequency. Thus, the roughness in 
the shapes of the emissivity spectra is another manifestation of 
systematic errors. 


B. Histograms of Emissivity Retrievals 

Further insight into the uncertainties can be obtained from 
histograms of the retrievals- Fig, 3 shows the histograms of 
horizontally polarized emissivities for both sites. Overall, the 
Desert site exhibits a gradual increase in H-pol emissivity with 
frequency (left) from all the retrievals, while over Amazon 2 
the emissivities are largely confined in the range of 0.9 to 
1 .0 for all the frequencies. Consistent with Fig. 2, there are 
considerable differences among the mean value of each of the 
retrievals over either site — a strong indicator of the existence 
of systematic errors. In addition, for each retrieval dataset, 
there is a range in spread around its mean emissivity value, 
wilh the shape of the histogram indicating the distribution of 
the random errors. 

Over Amazon2, all the emissivity histograms are single 
mode, while over Desert, some of the retrievals, such as I ’M I 
and AMSR-E (MIRS), exhibit dual modality at some of the 
frequencies. We speculate this might be related to the strong 
diurnal cycle in the variation of the surface temperature and 
the microwave penetration depth, which makes it tricky to 
represent the effective emission characteristics in the retrieval 
process [9], [18], [32], [34]. 

C. Uncertainties in MPDI 

Because MPDI can largely cancel the effect of errors in 
atmospheric contribution and surface temperature, we also 
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Fig. 3. Intercom paris on of the histograms of the H-pol microwave emissivity retrievals at two LSWG sites Desert (left) and Amazon 2 (right), over a one-year 
period from July 1, 2006 to June 30, 2007, for SSM/I (F13, F14, and F15), TMI, and AMSR-R, for various frequencies, (a) 10.65G. (b) 18.7/19.35G. (c) 
23. SG. (d)36.5/37.GG. (e) S5.5/89.GG, (D I0.65G. (g) 1 8.7/1 9.35G. (h) 23.8G. (i) 36.5/37.0G. (j) 85.5/89.0G. 


examined the uncertainties in MPDI computed from the emis- 
sivity retrievals. This will reveal how much of the error is 
common in both V- aud H-pol channels, and how much is not. 
if a significant portion of the error in an emissivity retrieval 
is common to both channels, MPDi will show much lower 
systematic and random error amplitudes than the emissivity 
values alone. Indeed, as shown in Fig. 4{a), MPDI values 
over Desert show much better agreement among the retrievals* 
except AMSR-E (M1RS) + This suggests most of the system- 
atic errors in the emissivity retrievals are common to both 
channels, in addition, the variance becomes much smaller, 
indicating that both channels have the same random error as 
well most of the times. Similar conclusions can be drawn for 
Amazon2 [Fig. 4(c)] for most of the frequencies, except for 
85.5/89.0 GHz, which shows fairly large systematic and ran- 
dom errors. Obviously, the errors at ihe highest frequencies are 
larger and less co- varying between the vertical and horizontal 
polarizations. 

For comparison, we also studied the MPDi computed 
from AMSR-E TOA Tb values over these two sites 
| Fig. 4(h) and -(d) J- Interestingly, the spectral shapes 
between emissivity-based (Fig. 4, left) and Tb-based MPDI 
(Fig. 4, right) are strikingly similar. Since Tb-hased MPI>1 is 
the lower bound of emissivity-hased MPDI, any values in the 
former (Fig. 4, right) higher than those in the latter indicate 


systematic errors in the emissivity retrievals, such as the FI 5 
retrievals at higher frequencies over Amazon2 [Fig. 4(c)]. In 
addition, there is an elbow in Tb-based MPDI at 23.8 GHz 
over either site, due to the strong water vapor attenuation in 
the atmosphere. But such a depression should not be presenf in 
emissivity-based MPDI, had the atmospheric effect been com- 
pletely removed in the retrieval process. But its very existence 
[Fig, 4(a) and (c)] suggests otherwiscn This is understandable 
because there is no way to recover the polarization difference 
from the TOA Tb if such a signal is strongly dissipated through 
the atmosphere. Therefore, satellite -based direct retrieval of 
land-surface emissivity in this frequency range (21-24 GHz) 
will remain a challenge, and a feasible solution is interpolation 
from its more transparent neighboring frequencies, as done 
by the tool to estimate land-surface emissivities at microwave 
frequencies (TELSEM) [35] . 

D. Error Diagnosis 

To further understand the causes of the errors, we inspected 
the daily emissivity spectra and their variations over our 
l -year period, for the various sensors for both sites. As 
an illustrative example, Fig. 5 shows a collage of daily 
emissivity retrievals from SSM/1 on DMSP FI 3 over the 
Amazon2 site, for both polarizations and both morning 
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Fig. 4. Intercomparison of emissivity -based MPDI values (left) at two LSWG sites Desert (top) and Amazon2 (bottom), over a one-year period from July 
I, 2006 to June 30, 2007, for SSM/T (FI 3, FI 4, and FI 5), TMT, and AMSR-E. For comparison, AMSR-E Tb-based MPDI values over the same two sites are 
also shown (right), (a) Desert, e mis sivity- based MPDI. (b) Desert, Tb-based MPDI. (c) Amazon2 J emissivity- based MPDI. (d) Amazon2 t Tb-based MPDI. 


(A.M.; descending) and afternoon (P.M.; ascending) passes. 
The emissivities for both polarizations are of similar values, 
as expected over such a site. Most of the emissivity spectra are 
confined within the range of 0.9-1 .0, but there is considerable 
variation, driven largely by the seasonably. This variation 
is more likely introduced by the errors in bolh land-surface 
temperatures and atmospheric profiles used in the retrievals. 
The impact of water vapor attenuation at 22.2 GHz for the 
vertical polarization, in the shape of an elbow, from either 
A.M. or P.M. passes, is obvious J Fig. 5(a) and (c)|. Such an 
elbow does not manifest itself in the horizontal polarization 
due to the lack of the 22.2 GHz |Fig. 5(b) and (d)J. 

Fig. 5 also reveals that the dominant source of random errors 
at the highest frequency (85,5 GHz) is likely the contamination 
from cloudy or rainy skies. As one can observe, during the 
morning passes, there were very few oulliers at 85,5 GHz, But 
the afternoon passes saw a significant number of oulliers with 
drastically reduced emissivity values in both polarizations at 
this frequency [Fig. 5(c) and (d)j + This coincides with ihe rainy 
time of ihe Amazon precipitation diurnal cycle [36J, and the 
outliers indicate strong scattering from ice particles aloft. This 
suggests that the cloud screen step in the retrieval processes 
missed a number of cloudy and rainy conditions, resulting in 
enhanced random errors. 


E. Analysis of Error Sources 

There are numerous sources contributing to both the sys- 
tematic errors and the random errors. They are identified and 
their relative importance is discussed as follows. 

1) Differences in each of the sensors" configurations and 
characteristics. These include differences in the fre- 
quency, incidence angle, footprint size, overpass time, 
and errors in brightness temperature measurements. 
However, over land surfaces, the slight differences in 
frequencies and incidence angles among the imagers 
will not cause appreciable discrepancies in the retrieved 
emissivities, due to their weak dependencies on these 
two variables [23]. Also under our assumption that 
the radiometric properties over either the Desert or the 
Amazon site are sufficiently uniform and constant, the 
differences in footprint size and overpass Lime should 
not be considerable. The very small differences among 
the MPDI values (except M1RS; Fig. 4) support this 
assumption. In fact, even over many other types of 
surfaces, the emissivities from these different platforms, 
when retrieved with the same method and same ancillary 
datasets, are very close [37 1. In addition, many existing 
calibration/validation (Cal/Val) studies 1 1 9 hi 27 1 suggest 


Emissfvity Emissivity 


This article has bovn accepted Ibr inclusion in a Future issue of this journal. Content is Imal as presented, with Ihe exception of pagination. 


K 


fhiKK TRANSACTIONS ON CibUSCIKNCK AINU KKMUIfc SENSING 




Tig. 5, Daily microwave emissivity retrievals over Aniazon2 (bottom), over a one-year period from July 1, 2006 to June 30, 2007, for SSM/I F13. Both 
vertical (left) aad horizontal (right) polarizations, as well as descending (A,M. t top) and ascending (p. M . , bottom) passes, are shown. Each daily retrieval is 
designated hy a colored line. (a)-(d) Frequency (GHz). 


that the systematic differences between these sensors 5 
brightness temperature measurements are in the order 
of 1-3 K, much less than the error amplitude shown 
here. Therefore, we can infer that these factors do not 
dominate Lhe errors revealed in our analysis. 

2) Differences in retrieval methods. Two classes of methods 
are used in the retrievals; one is the variational method 
employed by MIRS, and the other is the direct solution 
of (1} and (2), used by the rest of the datasets. However, 
the retrieved variables from MIRS, including emissivity, 
will automatically satisfy (1) and (2). Therefore, if all 
the ancillary data are the same, theoretically both classes 
should produce identical results. A meaningful future 
effort would be testing these different retrieval methods 
with the same ancillary data, to verify Ibis conclusion. 


3) Differences and errors in ancillary data. Used to solve 
(1) and (2) for clear skies, these data include: rain/cloud 
masks, surface temperatures, and atmospheric temper- 
ature, and water vapor profiles. Yang and Wcng [15] 
studied the sensitivities of emissivities retrievals to 
these factors wilh specified errors, assuming clear skies. 
Our results indicate thal cioud/precipitation contamina- 
tion plays a significant role in practice (Fig. 5), even 
with best-effort rain/cloud screening. In addition, the 
considerable systematic differences seen here, even at 
lower frequencies, suggest that the discrepancies in 
the respective ancillary datasets are even larger than 
[15] tested, albeit in the current it is impossible to 
cleanly separate the contributions from these ancillary 
variables. Over the Desert, the impact of the error in 



